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Abstract 

The nucleon polarizabilities ue and (3m are studied in the context of the 
perturbative chiral quark model. We demonstrate that meson cloud effects 
are sufficient to explain the electric polarizability of nucleon. Contributions of 
excite quark states to the paramagnetic polarizability are dominant and cancel 
the diamagnetic polarizability arising from the chiral field. The obtained 
results are compared to data and other theoretical predictions. 
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I. INTRODUCTION 



The nucleoli polarizabilities a^v and (3^ with N = p or n are two of the fundamen- 
tal observables which encode the electromagnetic structure of the nucleon. They can 
be determined from the differential cross sections of Compton scattering off the nucleon 
7 + ^ 7 + iV (for a review see Refs. [1]- [5]). The idea to use Compton scattering to in- 
vestigate the internal structure of the nucleon was formulated more than 50 years ago [6,7]. 
The first detailed measurement of the nucleon polarizabilities was performed in 1960 [8] 
by analyzing elastic 7p scattering at 40-70 MeV. Compton scattering off the nucleon is 
constrained at very small energies by the Low-Gell-Mann-Goldberger low-energy theorem 
derived in Refs. [9,10]. It completely fixes the leading and the next-to-leading 0{u;) 
structure-independent contribution to the expansion of the scattering amplitude in powers of 
the photon frequency uu. The leading term is known as the Thomson term and is expressed 
in terms of electron charge and nucleon mass only. The second term is proportional to 
the square of the anomalous magnetic moment of the nucleon. The nucleon polarizabilities 
characterize the quadratic terms in u in the Compton scattering amplitude [11], of which 
the first structure-dependent terms come from Rayleigh scattering. Polarizabilities express 
the deformation of particles under the influence of quasi-static electric and magnetic fields. 

— * — * 

In particular, in the presence of an external electric E or magnetic B field the quark dis- 

— * — * 

tribution of the nucleon becomes distorted, leading to the induced electric d — AiraEE or 
magnetic jl — Att^mB dipole moments, where the coefficients aE and /9m ai'c the electric 
and magnetic polarizabilities of the nucleon. The interaction of the electric d, and magnetic 
jl dipole moments with the electromagnetic field leads to a change of the corresponding 
effective Hamiltonian. It describes the interaction of nucleons with electromagnetic fields 

— * — * 

generating quadratic terms in E and B: 

AH,fi^-—{aEE' + PMB'). (1) 

The nucleon polarizabilities can be extracted from the Compton scattering cross section off 
the nucleon by comparison with 

^r^[(aE + PM){l + cosef + (q;e-/?m)(1-cos^)'] + ... (2) 

where 9 is the scattering angle. It was proposed already in the 50s and later indicated in 
experimental and theoretical studies that the meson cloud effects give a significant contri- 
bution to the nucleon polarizabilities [5]. The sum aE + is constrained by the so-called 
Baldin sum rule [12] which establishes a connection to the integral over the unpolarized 
photo-absorption cross section. The current experimental data for nucleon polarizabilities 
are [5]: 

= 12.0 ± 0.6 X 10"^ fm^ , = 1.9 T 0.6 x 10"^ fm^ , 
al = 12.5 ± 1.7 X 10"^ fm^ , = 2.7 ^ 1.8 x 10"^ tm^ . 

Theoretical analyses of nucleon polarizabilities were done using different phenomenological 
approaches [13,14] (quark, soliton and Skyrme models) and methods of chiral perturbation 
theory [15-17]. 
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A nonrelativistic approach in second-order perturbation theory gives the well-known sum 
rules for cce = dig -|- Act^; and (3m — Pm'^'^ + P'm'j where for example 



a 



2 2. ET— TT , = -— (3) 







are the leading contribution and the retardation correction to the electric polarizability. The 
magnetic polarizability contains the para- and diamagnetic contributions with 

Here | n > and are eigenstate and energy of an excited state, d is the electric dipole 
moment, dz and /i^ are the -^-components of the electric and magnetic moments, respectively, 
Zjv and uin are the charge and the mass of nucleon, Qi, rrii, and pi are the charge, mass 
and internal coordinate of the constituent quark, respectively. The problems associated with 
the use of the nonrelativistic sum rules have been discussed in [2,5,16-18]. In particular, 
the proton electric polarizability is predicted to be larger than the neutron one, which is in 
conflict with the central values of the experiment with > [16]. It turns out that the 
inclusion of the A(1232)-resonance is essential in the calculation of Pm- It gives a large posi- 
tive contribution to the paramagnetic moment which compensates a relatively large negative 
diamagnetic polarizability [16,18]. Finally, the main problem is that the above nonrelativis- 
tic definition of the nucleon polarizabilities is not identical to the physical observables 
and Pm because of missing relativistic corrections [2,5]. Therefore, an accurate analysis of 
nucleon polarizabilities should be done within a truly relativistic approach. 

Assuming P— , C- and T-invariance, the general amplitude for Compton scattering can 
be written in terms of six real structure dependent functions Af{u, 9) with i = 1 . . . 6, where 
uj = uo' denotes the photon energy in the center-of-mass (CM) frame. The latter is specified 
by the following choice of momenta of initial [p) and final {p') nucleons and incoming [k) 
and out-coming {k') photon: 

p = {E, -k) , p' = {E, -k') , k={u;,k), k' = {oj, k') . (5) 

For the polarization vectors of the real photon field it is convenient to choose the Coulomb 
gauge which is characterized by the conditions: 

e = (0,e), e*' = (0,e*'), A^-e = 0, k' ■ T' = Q . (6) 

The T-^-matrix of Compton scattering is then written as [19] 

= ^ Af(a;,^)5,, S^^e*'-e, S2 ^ e*' ■ k e- k' , ... (7) 

i=l 

Here, e. A; (e*', k') are the polarization vector and the direction of the incoming (outgoing) 
photon with k = k/\k\. The scattering angle 6 is defined as cos^ = k ■ k' . In Eq. (7) only 
the first two terms are spin-independent and others are spin-dependent ones which contain 



3 



the nucleon spin operator a. The structure functions contain contributions from the 
pion-pole (anomalous) and from the remaining (regular) terms with 

Af{uj, e) = (a;, ^)-°-P°'- + Af{uj, Of'^, i = 1...6. (8) 

However, the vr^'-pole does not contribute to the spin-independent structure functions A^ 
and A2 . One can perform a low-energy expansion of the six nucleon structure functions 
Af{u, Oy^^ in powers of the photon energy cv. The exphcit expressions of the expansion for 
the case of a proton target of mass with anomalous magnetic moment can be found 
in Refs. [3,5]. The explicit expressions of the spin- independent functions A^(u},9y'^^ and 
A2 {uj, 9y^^, which are expressed in terms of as and /3m, are [15] 

A^{uj,ey''^^- — Ql + 47r(a^ + cos9/3Z)uj^ 

ITlN 

+ — (a^ + PZ) (1 + cos 0)u;' + O L\ -\-) , 

A^{u;,9y^^ ^^Q%u; - AnP^u;' - —{a^ + PZ)^' + oL\ ^) , (9) 
mjf ruN V 111,% J 

where Qtv is the nucleon charge. 

In this work, a perturbative chiral quark model (PCQM) [20-22] is employed to study the 
spin-independent polarizabilities of nucleon, aN and Pn- In the PCQM baryons are described 
by three relativistic valence quarks confined in a static potential, which are supplemented 
by a cloud of pseudoscalar Goldstone bosons, as required by chiral symmetry. This model 
has already been successfully applied to the charge and magnetic form factors of baryons, 
sigma terms, ground state masses of baryons, the electromagnetic N ^ A transition, and 
other baryon properties [20-22]. A detailed discussion of the approach and a comparison to 
similar chiral quark models can be found in Ref . [20] . To calculate the nucleon polarizabihties 
we apply the heavy nucleon mass limit ^ 00 where only the cu^ terms survive in the 
expansion of the regular amplitudes Af^-^ 2{'^ 1 ^Y^^ ■ From an evaluation of these terms we 
extract the nucleon polarizabilities. 

The paper is organized as follows. In Sec. II we present an introduction to the PCQM. 
In Sec. Ill we calculate the nucleon polarizabilities in the PCQM. We discuss the numerical 
results and compare them to the predictions of other theoretical approaches. In Sec. IV we 
give our conclusions. 



II. THE PERTURBATIVE CHIRAL QUARK MODEL 

The perturbative chiral quark model is based on an effective chiral Lagrangian describing 
baryons by a core of the three valence quarks, moving in a central Dirac field with Kff('") — 
S{r) + 7°K(r), where r =\ x \. In order to respect chiral symmetry, a cloud of Goldstone 
bosons (tt, K and rj) is included as small fluctuations around the three-quark core. Here 
we restrict to the SU(2) version of our approach and take into account only the pion cloud 
fluctuations. The total model Lagrangian is 

-^eflf — Anv + i^xSB, (10) 
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where is the chiral-invariant Lagrangian and C^sb is the chiral symmetry breaking part. 
Up to second order in the pion field tx — Tr^Ti fiuctuations, >Cinv and jC-^^sb are 

Anv(x) = V'(x) [t - S{r) - 7V(r)]V^(a:) + ^id^Hx)]' - ^{x)S{r)i^'^il:{x) (11) 



^xSb{x) = -7P{x)M^{x) - ^Tr[e{x)M] . (12) 



Eq. (12) contains the mass contributions both for the quarks and pions, which exphcitly 
break chiral symmetry. Lagrangian (11) contains the basic parameters: F = 88 MeV is the 
pion decay constant in the chiral limit, — diag{m, rh} is the mass matrix of current quarks 
(we restrict to the isospin symmetry limit = = m) and B — — < 0\uu\0 > /F^ = 
— < 0|(irf|0 > /F^ is the quark condensate constant. We rely on the standard picture of 
chiral symmetry breaking [23] and for the masses of pseudoscalar mesons we use the leading 
term in their chiral expansion (i.e. linear in the current quark mass). By construction, 
our effective chiral Lagrangian is consistent with the known low-energy theorems (Gell- 
Mann-Okubo and Gell-Mann-Oakes-Renner relations, partial conservation of axial current 
(PCAC), Feynman-Hellmann relation between pion-nucleon cr-term and derivative of the 
nucleon mass, etc.). 

We expand the quark field ip in the basis of potential eigenstates as 

a 

where the sets of quark {ua] and antiquark {vp) wave functions in orbits a and P are 
solutions of the Dirac equation with the static potential Vef[{r). The expansion coefficients 
ha and rf^ are the corresponding single quark annihilation and antiquark creation operators. 
In our calculation of matrix elements, we project quark diagrams on the respective baryon 
states. The baryon states are conventionally set up by the product of the SU(6) spin-flavor 
and SU(3)c color wave functions, where the nonrelativistic single quark spin wave function 
is simply replaced by the relativistic solution Ua{x) of the Dirac equation 

' -i-i^^ ■ V + 7°^(r) + V{r) - 8^] u^ix) = 0, (13) 

where Sa is the single-quark energy of state a. 

For the description of baryon properties, we use the effective potential Ves{r) with a 
quadratic radial dependence [20,21]: 

S{r) = Ml + Clr^ V{r) = M2 + car' (14) 

with the particular choice 

M 1 - p 1 + _ p 
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Here, So is the single-quark ground-state energy; R and p are parameters related to the 
ground-state quark wave function uq: 



uo{x; i) = No exp 



X 



2 



where A^o = [7r^/^i?^(l+3p^/2)]~-'^/^ is a normalization constant; Xs: Xf: Xc are the spin, flavor 
and color quark wave functions, respectively. The index "i" stands for the i-th quark. The 
constant part of the scalar potential Mi can be interpreted as the constituent mass of the 
quark, which is simply the displacement of the current quark mass due to the potential S{r). 
The parameter p is related to the axial charge of the nucleon calculated in zeroth-order 
(or 3q-core) approximation: 



Therefore, p can be replaced by qa using the matching condition (17). The parameter R is 
related to the charge radius of the proton in the zeroth-order approximation as 

<rl>io= j d^x4{x)x^Uo{x) = ^ (18) 

In our calculations we use the value §'^=1.25. Therefore, we have only one free parameter 
in our model, that is or < r% >lo- the numerical studies, R is varied in the region 
from 0.55 fm to 0.65 fm, which corresponds to a change of < r|; >lo from 0.5 to 0.7 fm^. 
Note that for the given form of the effective potential (14) the Dirac equation (13) can be 
solved analytically [for the ground state see Eq.(16), for excited states see Ref. [22]]. In the 
Appendix we give details of the solutions to the Dirac equation for any excited state. 
The expectation value of an operator A is then set up as: 



CO -n I- I- 

MY.- d'x, ... dVr[£,(xi) . ..Cj{xr,)A]\<t>o >f , (19) 

n=l •' 



where the state vector |0o > corresponds to the unperturbed three-quark state (3g-core). 
Superscript "i?" in the equation indicates that the matrix elements have to be projected onto 
the respective baryon states, whereas subscript "c" refers to contributions from connected 
graphs only. Ci{x) of Eq. (19) refers to the quark-meson interaction Lagrangian: 

C:{x) = -V5(x)if ^5(r)V(x) + i^{xf^S{T)^l:{x) . (20) 

For the evaluation of Eq. (19) we apply Wick's theorem with the appropriate propagators 
for the quarks and pions. 

For the quark field we use a Feynman propagator for a fermion in a binding potential 
with 

iG^{x,y) = < (po\T{^|J{x)^|J{y)}\(f)o > 

= 0{xo - ?/o) E«a(^)^a(y)e-^^«(^°-^°) - 9{yo - xo)Y.vp{x)vM^'''^'''''°^ ■ (21) 
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In Refs. [20,21] we restricted the expansion of the quark propagator to its ground state with: 

iG^{x, y) ^ iGo{x, y) = uo{x) uo{y) e-*^o(-o-^o) e{xo - yo). (22) 

In Ref. [22] we investigated in addition the impact of excited states in the loops on the 
electromagnetic properties of baryons, the N-A transition, the nucleon axial vector form 
factors, and meson-baryon sigma-terms. Here, we also include the excited quark states in 
the quark propagator of Eq. (21) to analyze their contribution to the nucleon polarizabilities. 
In our numerical calculation we include the following set of excited quark states: the first 
p-statcs (lpi/2 and lp3/2 in the non-relativistic notation) and the second excited states 
{1(1^/2, 1(^5/2 and 2S1/2), i-e. we restrict to the low-lying excited states with energies smaller 
than the typical scale of A = 1 GeV. The justification for such an approximation can be 
found in our previous papers [22]. 

For the pions we adopt the free Feynman propagator with 

iA,j{x - y) =< 0|T{7r,(x)7r,(y)}|0 >= 5,, j ^ e-^(^-^) A.(A;) , (23) 

where A^{k) = [M^ - P-iO^]-^ is the pion propagator in the momentum space. Introduc- 
tion of the electromagnetic field ^4^ to the effective Lagrangian is accomplished by minimal 
substitution: 

Dfj^ijj = + ieQAfj^iJj, ^^tt^ D^iii = d^-Ki ee^ijA^/Kj (24) 

where Q — diag{2/3, — 1/3} is the SU{2) quark charge matrix and Sijk are the totally 
antisymmetric structure constants of SU{2) (Levi-Civita tensor). 
For the photon field A^ we also include the usual kinetic term 

Cph^-\F^.F>'- with F^,^d^A,-d,A^. (25) 
The electromagnetic interaction of the system is 

Cem^-eA^i^-tl'Q^P-e A^ e^j 7^8^" ttj + -A^A^ (tt^ - tt^) . (26) 



III. NUCLEON POLARIZABILITIES IN THE PCQM 



In this Section we apply the PCQM to calculate the spin-independent polarizabilities as 



and (3m contained in the real Compton scattering amplitudes A^ 



N, 



In the heavy nucleon mass limit with rriN 
amplitudes are sufficiently simplified as: 



u,ey''^ and A^{u, 

Ni 



\reg 



00 the structures of A^{uj, 6)'^^ and Ag{uj, 6^"^ 



A^{u,ey''^ 
A^{u;,ey^^ 



miv— >oo 



(27) 



miv— >oo 
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This limit it well-justified in the PCQM and corresponds to the static three-quark bag with 
projection onto a nucleon state with finite three-momentum. Let us stress again that we do 
not intend to calculate the full Compton scattering amplitude j + N —>■ ■y + N consistent 
with the relativistic theory. We just restrict to the evaluation of the spin-independent 
polarizabilities ajv and /3n- In the calculation of ajv and we proceed as follows. First 
we calculate the quark (one- and two-body) operators describing Compton scattering off a 
quark: 7 + 5-^7 + 5. Then we project the set of these operators onto the nucleon states 
and extract and /3jv from matching the result to Eq. (27). 

The Feynman diagrams contributing to the process 7 + g — > 7 + g (Compton scattering 
on the quark level) are shown in Figs. 1-8: tree-level diagrams (Fig.l), one-body meson- 
cloud diagrams (Figs. 2, 3, 4, 5 and 6) and two-body meson-cloud diagrams (Figs. 7 and 
8). Note, that the diagrams in Figs. 4, 5, and 7 do not contribute to the real part of the 
spin-independent polarizabilities of the nucleon aE and Pm- They only contribute to the 
spin-dependent polarizabilities. The intermediate ground state quark propagators contribute 
only in the diagrams of Figs. 2 and 8 (both to ccg and to Pm)- In the case of the diagrams 
of Fig. 6 the intermediate ground quarks state contribute to the magnetic polarizabilities 
only. Excited intermediate quark states contribute to the tree- level (Fig.l) and meson-cloud 
diagrams (Figs. 2 and 6). 

For illustration we present the expressions for some typical diagrams contributing to the 
matrix. The contribution of the tree-level diagrams (Figs. la and lb) is: 

rpN */ Tmn jmn jmn _i_ jmn 

^1 -"l ~-'la"r-'lb' 

rrnn _ 2 ^ Sr YLlih^l}^^^ .Tmn _ o ' ^ Vrj-t, a) Vr\-k ' ; a) 

where {m,n) — 1,2,3, and AEq, is the energy difference between the intermediate quark 
state (labeled by a) and the initial (or final) quark ground state. The vertex function of the 
photon-quark coupling in Eq. (28) is 

V;'\k; a)= J d^xuoix; th"VtUa{x; ^)e^*^'*, (29) 

where UQ{x;i) and Ua{x;i) are the wave functions of ground and excited quark states; Qi 
is the charge number of the i-th quark. The amplitude /{"" is proportional to J2i=i Qj, 
resulting in 1 for proton and 2/3 for neutron. If only the ground state is considered, the 
three-component of the vertex function is 

Vi(k; 0) = iN^pRinRY^^e-^ Qi [a(i) x k], (30) 

of which the quantities R, p, and Nq are defined in Sec. II. Tree-level diagrams constrained to 
include only intermediate ground state quarks do not contribute to the nucleon polarizabil- 
ities. This is due to the fact that A.Ea — and, therefore, the expansion of T{"" in powers 
of cu does not contain terms. Inclusion of the intermediate excited quark states in the 
tree-level diagrams leads to a contribution to the nucleon polarizabilities due to AE^ 7^ 0. 
Moreover, only /3m is affected by the tree-level diagrams, because they have the following 
spin structure 
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a{i) ■ ik' X e *') ■ {k x e) — a;^(cos 9 e ■ e*' — e*' ■ ke • k') + spin — dependent part . 

As already mentioned before, only the meson-cloud diagrams of Figs. 2, 3, 6, and 8 
contribute to the spin- independent nucleon polarizabilities q;^ and (5m- In other chiral 
approaches, the one-body diagrams of Figs. 2 and 3 give rise to the dominant terms to 
ue and (3m- In addition, in our approach we take into account two-body quark forces. 
Therefore, ue and (^m receive a contribution from the diagrams of Fig.8 which arc the two- 
body analogue of the one-body diagrams of Fig.2. The diagrams of Fig. 6 evaluated for the 
nucleon polarizabilities are suppressed. For example, the explicit expression of the diagrams 
in Fig.2 to the matrix is given by 

rpN */ jmn jmn jmn , jmn , jmn ('^^\ 

-'-2 —^m^2 ^n-i ^2 —-'2a ^2h ^2c i I'^-'-J 

where 

e2 3 ^4^ 



jmn I jmn ^ \ ^ \ ^ / 

-'2a + -'26 — p2 2^ 2^ J 



i=l a 

X " 



{2t:)H 



P^{q-k';z)A^{q-k')A^{q)A^{q-k)Pl{q-k;z)^^'^ ^''^'"^"'^ 



qo + AEc, -uj-i0+ 

+ Pa{q + k i)A^{q + k')A^{q)A^{q + k)Pl{q + k'; i) ' 



jmn Xi^iT' ^ \ ^ \ ^ / 



qo + AEa+u;- iO+. 

(32) 

e2 3 1 



{2tt)H q^ + ^ + AEa- tO+ 



„ Ak \ ^ / Ak\ ^ ( Ak\^.(^ Ak 

— * 

with q being the pion momentum and Ak — k' — k. Here, Pa{k; i) is the quark-pion vertex 
function 

Pa{k; i)^ j d^xuoix; i)^5S{r)ua{x; iy^'^ (33) 

where S{r) is the scalar part of the confinement potential. In particular, when a = 0, we 
have 

Po{k- i) = QA F.iviv(^') ■ k, (34) 

^ ,r2N ( k'^R^\ r k^R^ ( 5 \ 
W*^)=exp(-— ) + 

where FTjNN{k'^) is the ttNN form factor normalized to unity at zero recoil, fc^ = 0. 

Our results for the proton and neutron spin-independent polarizabilities and Pm 
are contained in Figs. 9- 18. In Figs. 9-12 we show the sensitivity of the proton (Figs. 9 and 
11) and neutron polarizabilities (Figs. 10 and 12) to a variation in the model parameter 
< r|; >Lo- I^ Figs. 13-18 we demonstrate the contribution of the individual excited states in 
the tree-level (Figs. 13 and 14) and in the meson-cloud (Figs. 15-18) diagrams to (3^^ and 
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The contribution of the tree-level diagrams to /3m grows linearly when the size parameter 
< r|. >LQ increases. Results for the meson-cloud diagrams are not changed so much as 
function of < >lo (when we restrict to the inclusion of the ground state quark propagators 
only) and grows (decreases) in (/3m) when < r"^ >^q increases. 

Now we discuss the numerical results for the nucleon polarizabilities in the case of the 
central value < r"^ >lo~ 0-6 f^n^ obtained in previous phenomenology [20-22]. In Table I 
we compare our results to the data and to other theories (ChPT [15], ChPT -|- dispersion 
relation (DR) constraints [17]). In Table II we present the individual contributions of the 
different diagrams to the nucleon polarizabihties, while also hsting the separate results for 
the ground and excited quark states in the loops. It can be seen that the tree-level diagrams 
and the meson cloud diagrams of Fig. 6 do not contribute to the electric polarizabilities. The 
contribution of the tree-level diagrams to the magnetic polarizabilities is solely due to the 
inclusion of the excited quark states. As already indicated, the excited state contributions 
in the tree-level diagrams are for the case of the neutron 2/3 of the result of the proton. 
The total predicted value for /3^ is larger than that of /3". Finally, we found that only very 
small differences appear in the contribution of the six diagrams of Fig. 6 to /?m, which changes 
slightly from 0.085 for the proton to 0.089 for the neutron. In our model calculations there is 
no remarkable difference between the proton and neutron cases for the pion-cloud diagrams. 
Thus, our results show that ctg = a^. Moreover, the large difference between proton and 
neutron arises for Pm, where the tree-level diagrams of intermediate excited quark states 
differ strongly for both cases. 

IV. CONCLUSIONS 

In this work we apply the perturbative chiral quark model to the description of the 
nucleon spin-independent polarizabilities as and Pm- It has been previously verified that 
the model is successful in the explanation of many aspects of nucleon properties [14-18], such 
as the magnetic moments, the axial vector form factor, the N ^ A transition amplitude, 
the meson-nucleon sigma-term and nN nucleon scattering. Here, we find that the model can 
reasonably reproduce the electric polarizabilities of proton and neutron without the use of 
any additional free parameters. The magnetic polarizabilities, in particular for the proton, 
are overestimated when compared to the present data. This effect is known from experience 
of other theoretical calculations. 

In the perturbative chiral quark model it is found that without excited quark states in 
the loops already reasonable value for and a negative one for /?m are obtained. Then 
(3m behaves like a diamagnetic polarizability. This finding also appears in the dispersion 
relation analysis, in other chiral quark model calculations and in leading-order relativistic 
chiral perturbative theory (see discussion in Ref. [5]). When we include the excite states in 
the quark propagators we get a paramagnetic polarizability (3m arising from the tree-level 
diagrams. Our value for (3^ is overestimated in comparison with the data. This trend is also 
seen in heavy baryon chiral perturbative theory with the missing A(1232) contribution in 
terms of a "small scale expansion" . From a detailed analysis we find that the intermediate 
excited quark states are important to generate a positive value for (3m- In fact, the sizeable 
contribution of excited states has already been elaborated in our calculation of the A —> N 
transition amplitude and in the case of the nucleon axial charge [22]. Comparing each 
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contribution of the five intermediate excited states, we find that the efi^ects of and 
lpi/2 are more pronounced than the effects of the other three resonances. For the neutron 
case, the main contribution, to the paramagnetic polarizabihty j3M from the excited states 
of the Born terms, is proportional to Z)|=i Thus, the paramagnetic polarizabihty of the 
neutron is smaller than that of the proton target. For the observable a^, the one-pion loop 
contributions, both from the ground state and from the excited states, are not sensitive to 
isospin. Consequently, we have a% ~ 0% in the present calculation. 
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APPENDIX A. SOLUTIONS OF THE DIRAC EQUATION FOR THE 

EFFECTIVE POTENTIAL 



In this section we indicate the solutions to the Dirac equation with the effective potential 
Ves{r) — S{r) + 7°y(r). The scalar S{r) and time-hke vector V{r) parts are given by 



S{r) = Ml + cir^ 
V{r) = M2 + C2r\ 



with the particular choice 

Ml : ^ 



M2 = £0 - 



1 + 3p2 



Ci = C2 



(35) 



(36) 



2pR ' '"^ 2pR ' 2i?3- 

The quark wave function Ua{r) in state a with eigen-energy Sa satisfies the Dirac equation 

i-iaV + (5S{t) + V{t) - 8^]u^{t) = 0. (37) 
Solutions of the Dirac spinor Ua{r) can be written in the form [24] 



For the particular choice of the potential the radial functions g and / satisfy the form 



(38) 



9a{r) 



Rn 



n-1 



\Rl) 



(39) 



where for j' = Z + ^ 



R^ 



(40) 



and ion j — I — \ 



fair) = -Pa( — 



r2 



(41) 



The label a = {nljm) characterizes the state with principle quantum number n = 1, 2, 3, 
orbital angular momentum I, total angular momentum j — I ± ^ and projection m. Due 
to the quadratic nature of the potential the radial wave functions contain the associated 
Laguerre polynomials L^{x) with 



^n(^) = E(-iy 



(n + A;)! 



^=0 {n — m)\{k + m)\m\ 

The angular dependence (3^a('^) = yimj{r)) is defined by 

yimj{r)^ {im^ms\jm)Yi,n^{r)x^ 

mi,ms 



X" 



(42) 



(43) 
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where Yijni{r) is the usual spherical harmonic. Flavor and color parts of the Dirac spinor 

are represented hj Xf and Xc, respectively. 

The normalization constant is obtained from the condition 



(44) 



which results in 



(2n + 2/)! 
""(n + OKn- 1)! 



2-2(n-^.+l/2)^l/2^3 ^ ^-—^ {1 + pii2n + I - ^)} 



-1/2 



The two coefficients Ra and pa are of the form 

Ra = R{l + ASapRy^^\ 
/R ^3 



Pa = P[ 



a 

V R 



(45) 

(46) 
(47) 



and are related to the Gaussian parameters p and R of Eq. (16). The quantity ASa = ^a~^o 
is the difference between the energy of state a and the ground state. A£a depends on the 
quantum numbers n and I and is related to the parameters p and it! by 



_P_ 

R^ 



(4n + 21-1) 



(48) 
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FIGURES 



Fig.l: Tree-level diagrams contributing to Compton scattering on a quark. 
Figs. 2-6: One-body meson-cloud diagrams in Compton scattering. 
Figs. 7-8: Two-body meson-cloud diagrams in Compton scattering. 

Fig. 9: Proton electric polarizability a^: The dotted and the dashed curves are the results 
of the pion loops without and with inclusion of the excited states. The total value is shown 
by the solid line. The shaded region shows the experimental results including error bars [5] . 

Fig. 10: Neutron electric polarizabihty a^: Otherwise as in Fig.9. 

Fig. 11: Proton magnetic polarizabihty (3^: The dotted and dashed curves describe the 
results of the pion loops without and with inclusion of the excited states, respectively. The 
dotted-dashed line is the contribution from the tree-level diagrams (excited intermediate 
quark states). The total value is shown by the sohd hue. The shaded region shows the 
experimental results including error bars [5]. 

Fig. 12: Neutron magnetic polarizability P^t- Otherwise as in Fig.ll. 

Fig. 13: Contribution of the tree- level diagrams (Fig.l) [excited intermediate quark states] 
to f3^. The dotted, dashed, dotted-dashed, two-dottcd-dashed and two-dashed-dotted lines 
stand for the results of the individual excited states lpi/2, lp3/2, 1(^3/2! 1(^5/2 and 2si/2. 
Fig. 14: Contribution of the tree- level diagrams (Fig.l) [excited intermediate quark states] 
to /3^. Otherwise as in Fig.l3. 

Fig. 15: Contribution of the meson-cloud diagram (Fig.2) [excited intermediate quark states] 
to a^. The dotted, dashed, dotted-dashed, two-dotted-dashed and two-dashed-dotted lines 

stand for the results of the individual excited states lpi/2, lp3/2, lc^3/2, 1^^5/2 and 2si/2. 
Fig. 16: Contribution of the meson-cloud diagram (Fig.2) [excited intermediate quark states] 
to a^. Otherwise as in Fig. 15. 

Fig. 17: Contribution of the meson-cloud diagram (Fig.2) [excited intermediate quark states] 
to The dotted, dashed, dotted-dashed, two-dotted-dashed and two-dashed-dotted lines 
stand for the results of the individual excited states lpi/2, lp3/2) 10^3/2, lc^5/2 and 2si/2. 
Fig. 18: Contribution of the meson-cloud diagram (Fig.2) [excited intermediate quark states] 
to P'^. Otherwise as in Fig. 17. 
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TABLES 



TABLE I. Results for the nucleon polarizabilities (in units of 10 ^fm'^). 



Approach 




0M 






data [5] 


12.0 ± 0.6 


1.9 T 0.6 


12.5 ± 1.7 


2.7 T 1-8 




7.9 


-2.3 


11.0 


-2.0 


ChPT [15] 


10.5 ±2.0 


3.5 ±3.6 


13.6 ± 1.5 


7.8 ±3.6 




12.6 


1.26 


12.6 


1.26 


ChPT + DR [17] 


7.3 


-1.8 


9.8 


-0.9 


PCQM 


10.9 


5.1 


10.9 


1.15 



TABLE II. Contributions to and f3^ from the individual diagrams (in units of 10 *^fm^) 
evaluated for the central value < r|; >lo= 0-6 fm^. 





a|(GS) 


a^(ES) 


a^(GS) 


a^(ES) 






/3g(GS) 




Fig.l 

















11.85 





7.90 


Fig.2 


2.03 


2.37 


2.03 


2.37 


-0.65 


-4.07 


-0.65 


-4.07 


Fig.3 


4.26 





4.26 





-1.42 





-1.42 





Fig. 6 














0.085 


0.02 


0.089 


0.02 


Fig.8 


2.26 





2.26 





-0.71 





-0.71 





Total 


8.55 


2.37 


8.55 


2.37 


-2.70 


7.80 


-2.70 


3.85 
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